First results of the search for neutrinoless double beta decay with the NEMO 3 

detector 
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The NEMO 3 detector, which has been operating in the Frejus underground laboratory since 
February 2003, is devoted to the search for neutrinoless double beta decay {fSpOv). The half- lives 
of the two neutrino double beta decay {(3[32u) have been measured for ^""^Mo and *^Se. After 389 
effective days of data collection from February 2003 until September 2004 (Phase I), no evidence for 
neutrinoless double beta decay was found from ~7 kg of ^°°Mo and ~1 kg of *^Se. The corresponding 
limits are T^/2{PPQv) > 4.6 x 10^^ years for i°°Mo and Ti/iiPPOu) > 1.0 x 10^^ years for ^^Se (90% 
C.L.). Depending on the nuclear matrix element calculation, the limits for the effective Majorana 
neutrino mass are {m„) < 0.7 - 2.8 eV for ^°°Mo and {m„) < 1.7 - 4.9 eV for *^Se. 



PACS numbers: 23.40.-s; 14.60.Pq 
I. INTRODUCTION 



The positive results obtained in the last few years 
in neutrino oscillation experiments Q, |^ |M Q have 
demonstrated that neutrinos are massive particles and 
that lepton flavor is not conserved. In parallel, tritium 
/3-decay experiments |^ have established a very low 
limit on the electron neutrino mass of m^^ < 2.2 eV 
(95% CL). Grand Unified Theories can provide a natural 
framework for neutrino masses and lepton number 
violation. In particular the see-saw model Q which 
requires the existence of a Majorana neutrino, naturally 
explains the smallness of neutrino masses. The existence 
of Majorana neutrinos would also provide a natural 
framework for the leptogenesis mechanism which 
could explain the observed baryon-antibaryon asym- 
metry in the universe. The observation of neutrinoless 
double beta decay {(3/30i') would prove that neutrinos 
are Majorana particles and that global lepton number 
is not conserved. It would also constrain the mass 
spectrum and the absolute mass of the neutrinos. 



II. THE NEMO 3 DETECTOR 

The NEMO 3 detector 0, installed in the Frejus 
underground laboratory (LSM, France) is searching 
for (3(30^ decay by the direct detection of the two 
electrons with a combination of tracking and calorimeter 
information. The two main isotopes present inside the 
detector in the form of very thin foils (40-60 mg/cm'^) 
are i"°Mo (6914 g, Qp^ = SOSikeV) and ^^Se (932 g, 
Q/313 = 2295keV). On both sides of the sources, there 
is a gaseous tracking detector which consists of 6180 
open drift cells operating in the Geiger mode allowing 
three-dimensional track reconstruction. To minimize 
the multiple scattering, the gas is a mixture of 95% 
helium, 4% ethyl alcohol, 1% argon and 0.1% water. 
The wire chamber is surrounded by a calorimeter which 
consists of 1940 plastic scintillator blocks coupled to 
very low radioactive photomultipliers (PMTs). The 
energy resolution (FWHM) of the calorimeter is 14% at 
1 MeV for the scintillators equipped with a 5" PMTs 
on the external wall and 17% for the 3" PMTs on the 
internal wall. The resolution of the summed energy 
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FIG. 1: Transverse (left) and longitudinal (right) view of 
a reconstructed (313 event selected from the data with a two 
electron energy sum of 2812 keV. 



of the two electrons in the P^Qv decay is mainly a 
convolution of the energy resolution of the calorimeter 
and the fluctuation in the electron energy loss in the 
foil source which gives a non-gaussian tail. The FWHM 
of the expected two-electron energy spectrum of the 
l3(3Qv decay is 350 keV. Absolute energy calibrations 
are carried out every 40 days using ^"^Bi sources. A 
daily laser survey controls the gain stability of each 
PMT. A solenoid surrounding the detector produces a 
25 Gauss magnetic field in order to distinguish electrons 
from positrons. An external shield of 18 cm of low 
radioactivity iron, a water shield and a wood shield 
cover the detector to reduce external 7 and neutrons. 



III. MEASUREMENT OF f3l32u DECAYS 

A two-electron (2e^) event (see Fig. ^ candidate for 
a 13(3 decay is defined as follows: two tracks come from 
the same vertex on the source foil, each track must be 
associated with a fired scintillator, its curvature must 
correspond to a negative charge and the time-of-flight 
must correspond to the two electrons being emitted from 
the same source position. For each electron an energy 
threshold of 200 keV for ^'^OMo and 300 keV for is 
applied. Fig. [2Ia) and Fig. |2| show the two-electron en- 
ergy sum spectra after background subtraction obtained 
after 389 effective days of data collection with ^°"Mo 
and with ^^Se respectively. The angular distribution of 
the two electrons and the single energy spectrum are 
also presented in the case of ^'^°Mo in Fig. |21[b) and (c). 
All these spectra are in good agreement with the f3(32v 
simulations. The values of the measured half-lives are 
T^/2{l3l32v) = [7.11 ± Qm{stat) ± 0.54(sysi)] x lO^^ y 
for ^"^Mo (with a Single State Dominance decay) and 
[9.6 ± 0.3(stat) ± l.O(sysi)] x lO^^ y for ^SSe. These 
values are in agreement with, but have a higher precision 
than the previous measurements 10]. 



IV. STUDY OF THE BACKGROUND IN THE 
I3I3QU ENERGY WINDOW 

A complete study of the background in the /3/90z/ en- 
ergy window has been performed. The level of each back- 
ground component has been directly measured using dif- 
ferent analysis channels in the data. 

External backgrounds due to ^^''Bi and ^osxi contami- 
nants outside the source foils (mostly in the PMTs) have 
been measured by searching for Compton electrons emit- 
ted from the source foils by external 7. For ^"^Tl, a total 
activity of ~40 Bq has been measured and is in agree- 
ment with the previous HPGe measurements of samples 
of the PMT glass. For 2i4Bi, an activity of -300 Bq 
has been found, again in agreement with the HPGe mea- 
surements of PMTs and also the level of radon surround- 
ing the detector inside the shield. The expected number 
of (3l3Qv-\\ke events due to this background is negligible, 
< 10-3 counts.kg-i.y-i in the [2.8 - 3.2] MeV energy 
window where the l3j3Qv signal is expected. 

External neutrons and high energy 7 backgrounds have 
been measured by searching for crossing electron events 
above 4 MeV. This corresponds to a negligible expected 
level of background of '--^ 3 x 10"^ counts. kg-^.y"^ in the 
PPOv energy window. 

The level of ^os^pj impurities inside the sources has 
been measured by searching for internal (6^77) and 
(e-777) events. The measured activity is 80±20 /iBq/kg 
in molybdenum and 300 ±50 /xBq/kg in Selenium. It is in 
agreement with the previous HPGe measurements which 
gave an upper limit of 100 /^Bq/kg for Molybdenum and a 
positive measurement of 400 ± 100 /iBq/kg for Selenium. 
This corresponds to an expected level of background 
in the PjSQv energy window of ~ 0.1 counts. kg-^.y^^ 
for Molybdenum and ~ 0.3 counts. kg^^.y"^ for Sele- 
nium. The measurement of ^^"^Bi impurities inside the 
sources could not be achieved in this first period of data 
due to Radon contamination (see later). However the 
previous HPGe measurements gave an upper limit of 
350 /xBq/kg for Molybdenum and a positive measure- 
ment of 1.2 ±0.5 mBq/kg for Selenium, corresponding to 
a negligable expected level of backgound. 

The expected level of background due to the tail 
of the P(32iy distribution in the /SpOv energy win- 
dow is ~ 0.3 counts. kg^^.y^^ for Molybdenum and ^ 
0.02 counts. kg- ^.y-^ for Selenium. 

The dominant background in this first period of data 
was Radon gas inside the tracking chamber due to 
a low rate of diffusion of Radon from the laboratory 
(^^15 Bq/m'^) into the detector. Two independent 
measurements of the Radon level in the detector were 
carried out. The first used a high sensitivity Radon 
detector similar to the one developed by the Super- 
Kamiokande collaboration [l^. The second was done 
by searching for (e-, delayed-a) events in the NEMO 3 
data. Indeed the tracking detector allows the detection 
of the delayed tracks (up to 700 /is later) in order to 
tag delayed-a emitted by ^^"^Po in the Bi-Po process. 
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FIG. 2: (a) Energy sum spectrum of the two electrons, (b) angular distribution of the two electrons and (c) single energy 
spectrum of the electrons, after background subtraction from ^"^Mo with 7.369 kg. years exposure. The solid line corresponds 
to the expected spectrum from /3/32i^ simulations and the shaded histogram is the subtracted background computed by Monte- 
Carlo simulations. The signal contains_2I9,000 /3/3 events and the signal-to-background ratio is 40. 
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FIG. 3: Energy sum spectrum of the two electrons after 
background subtraction from *^Se with 0.993 kg. years expo- 
sure (same legend as Fig. The signal contains 2,750 /3/9 
events and the signal-to-background ratio is 4. 



Both measurements are in good agreement and indicate 
a level of radon inside the detector of 25 ± 5 mBq/m^. 
This radon contamination corresponds to an expected 
level of background in the (3P0i> energy window of 
'--^l count. kg^^.y"^. 



V. FIRST RESULTS ON THE LIMITSEARCH 
OF ppOu DECAY WITH ^"MO AND ^^^SE 

Fig. EJa) and (b) show the tail of the two-electron en- 
ergy sum spectrum in the (3(3Qv energy window for ^"^"Mo 
and for ®^Se respectively. The number of 2e~ events ob- 
served in the data is in agreement with the expected num- 
ber of events from (3f32v and the Radon simulations. For 
i'^°Mo, in the energy window [2.8 - 3.2] MeV, the ex- 
pected background is 8.1 ± 1.3 (error dominated by the 
uncertainty on the Radon activity) and 7 events have 
been observed. For *^Se, in the energy window [2.7 — 3.2] 
MeV, the expected background is 3.1 ± 0.6 and 5 events 



have been observed. In order to independently check the 
dominant Radon contribution above 2.8 MeV, the energy 
sum spectrum (Fig.^c) has been plotted for the two elec- 
trons emitted from the Copper and Tellurium foils where 
no background except radon is expected. The data are 
in agreement with the Radon simulations. 

In order to make the optimum use of all the informa- 
tion from the NEMO 3 detector, a maximum likelihood 
analysis TTl has been applied to the 2e~ event sample 
above 2 MeV using the three available variables: the en- 
ergy sum {Etot) of the two electrons, the energy of each 
electron {Emin is the minimum electron energy) and the 
angle between the two tracks {cosO). A three-dimensional 
probability distribution function, P^^ , can be written as: 

P3D ^ P{Etot)P(E„^,n/Etot)P{cOSe/Ermn) 

where P{Ejnin/ Etot) and P[cos9 / Emin) are two condi- 
tional probability distribution functions. The likelihood 
is defined as 

where k corresponds to one of the eight contributions: 
l3(3Qv, I3j32v, Radon, external and internal •^^"'Bi and 
^"^Tl, and neutrons. Here Xk is the ratio of the num- 
ber of 2e~ events due to the process k relative to the 
total number of observed events Ntot- Finally P^^ is 
built using simulated events of contribution k. 

With 389 effective days of data collection, limits 
at 90% C.L. obtained with the likelihood analysis 
are ri/2(/3/30j/) > 4.6 x 10^3 years for i°°Mo and 
1.0 X 10^'^ years for ®^Se. These limits are about 10 
times higher than the previous limits obtained with 
^•^^Mo and ®^Se 0, The corresponding upper 

limits for the effective Majorana neutrino mass range 
from 0.7 to 2.8 eV for i"OMo and 1.7 to 4.9 eV for ^"^Se 
depending on the nuclear matrix element calculation 
15, 16, 17, 18, 19, 20] . Results for each calculation are 
given in Table |l| For ^''^Mo, since an incorrect value 
of the phase-space factor has been used in reference 
[20l |. the value calculated in reference 0| has been 
used. The claim of a positive /?/30j/ signal observed 
with ^^Ge |23| gives an allowed effective mass range 
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FIG. 4: Spectra of the energy sum of the two electrons in the PPOv energy window after 389 effective days of data collection 
from February 2003 until September 2004 (Phase I): (a) with 6.914 kg of ^°°Mo; (b) with 0.932 kg of ^^Se; (c) with Copper 
and Tellurium foils. The shaded histograms are the expected backrounds computed by Monte-Carlo simulations: dark (blue) 
is the /3/32i/ contribution and light (green) is the Radon contribution. The solid line corresponds to the expected PfSQu signal 
if ri/2(/3/30i/) = 5 X 10^^ y. 



TABLE I: Limits (in eV) on the effective neutrino mass (m^) 
obtained from different theoretical calculations of nuclear ma- 
trix elements with T^/iiPPQu) > 4.6 x lO^^y for "'"Mo and 
Ti/2(/3/30!^) > 1.0 X lO^V for '^^Se. 

Nuclear matrix elements ^°"Mo ®^Se 

Shell model Caurier 1996 fl5] <T9 

QRPA Rodin 2005 [16j < 2.7 - 2.8 < 4.1 - 4.5 

QRPA Simkovic 1999 [17] < 1.0 < 3.3 

QRPA Suhonen 2003 fl8, < 1.1 < 2.8 - 4.2 

QRPA Stoica 2001 [20] < 0.7 - 1.1 < 1.7 - 3.7 



0.1 - 0.9 eV. Our limit obtained with ^°°Mo slightly 
overlaps this range. In the hypothesis of gluino or 
neutralino exchange, and using the nuclear matrix 
elements calculated in reference |22 |. limits obtained on 
the trilinear R-parity-violating supersymmetric coupling 
are A'^i < 1-6 x IQ-* for io°Mo and 3.0 x 10"'' for ^^Se. 
In the hypothesis of a right-handed weak current, the 
limits are Ti/2{l3(3Qv) > 1.7 x 10^^ years at 90% C.L. for 
i00]y[Q Q 7 X 10^^ years for *^Se, corresponding to an 
upper limit on the coupling constant of A < 2.5 x 10^^ 
for lo^Mo and 3.8 x 10"^ for ^^Se using the nuclear 
calculations from references |0, . 



measured for ®^Se and ^"^Mo with very high statistics 
and better precision than the previous measurements. 
The two-electron energy sum spectrum, the single en- 
ergy spectrum and the angular distibution are all in good 
agreement with the j3[32v simulations. All components 
of the background in the PfSOv energy window have been 
measured directly using different analysis channels in the 
data. After 389 effective days of data collection, no evi- 
dence for PPOv decay has been found in ^"^Mo or ^^Se. 
The limits at the 90% C.L. are Ti/2{PP0v) > 4.6 x lO^^ y 
for ^"^Mo and 1.0 x 10^^ y for ®^Se. For this first running 
period (Phase I) presented here. Radon was the domi- 
nant background at a level of about 3 times higher than 
the f3l32i/ background for "'^"^Mo. It has now been sig- 
nificantly reduced by a factor '^10 by a radon-tight tent 
enclosing the detector and a radon-trap facility in op- 
eration since December 2004 which has started a sec- 
ond running period (Phase II). After five years of data 
collection, the expected sensitivity at 90% C.L will be 
Ti/2iPf30iy) > 2 X 10^4 y for ^''OMo and 8 x 10^3 y for 
^^Se, corresponding to (?7i^) < 0.3 - 1.3 eV for ^"''Mo 
and (m^) < 0.6 - 1.7 eV for ^^Se. 



VI. CONCLUSIONS 

A portion of this work was supported by grants from 
In conclusion, the NEMO 3 detector has been running INTAS (03051-3431, 03-55-1689) and NATO (PST CLG 
reliably since February 2003. The /3/32i/ decay has been 980022) 
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